Effects of Leishmania major infection on the gut microbiome of resistant and susceptible mice.

Cutaneous leishmaniasis, a parasitic disease caused by Leishmania major, is a widely frequent form in humans. To explore the importance of the host gut microbiota and to investigate its changes during L. major infection, two different groups of mouse models were assessed. The microbiome of two parts of the host gut-ileum and colon-from infected and non-infected mice were characterised by sequencing of 16S rDNA using an Ion Torrent PGM platform. Microbiome analysis was performed to reveal changes related to the susceptibility and the genetics of mice strains in two different gut compartments and to compare the results between infected and non-infected mice. The results showed that Leishmania infection affects mainly the ileum microbiota, whereas the colon bacterial community was more stable. Different biomarkers were determined in the gut microbiota of infected resistant mice and infected susceptible mice using LEfSe analysis. Lactobacillaceae was associated with resistance in the colon microbiota of all resistant mice strains infected with L. major. Genes related to xenobiotic biodegradation and metabolism and amino acid metabolism were primarily enriched in the small intestine microbiome of resistant strains, while genes associated with carbohydrate metabolism and glycan biosynthesis and metabolism were most abundant in the gut microbiome of the infected susceptible mice. These results should improve our understanding of host-parasite interaction and provide important insights into the effect of leishmaniasis on the gut microbiota. Also, this study highlights the role of host genetic variation in shaping the diversity and composition of the gut microbiome. KEY POINTS: • Leishmaniasis may affect mainly the ileum microbiota while colon microbiota was more stable. • Biomarkers related with resistance or susceptibility were determined in the gut microbiota of mice. • Several pathways were predicted to be upregulated in the gut microbiota of resistant or susceptible mice.

Thiothrix and Sulfurovum genera dominate bacterial mats in Slovak cold sulfur springs.

Microbiota of sulfur-rich environments has been extensively studied due to the biotechnological potential of sulfur bacteria, or as a model of ancient life. Cold terrestrial sulfur springs are less studied compared to sulfur-oxidizing microbiota of hydrothermal vents, volcanic environments, or soda lakes. Despite that, several studies suggested that sulfur springs harbor diverse microbial communities because of the unique geochemical conditions of upwelling waters. In this study, the microbiota of five terrestrial sulfur springs was examined using a 16 S rRNA gene sequencing. The clear dominance of the Proteobacteria and Campylobacterota phyla of cold sulfur springs microbiota was observed. Contrary to that, the microbiota of the hot sulfur spring was dominated by the Aquificota and Firmicutes phylum respectively. Sulfur-oxidizing genera constituted a dominant part of the microbial populations with the Thiothrix and Sulfurovum genera identified as the core microbiota of cold sulfur terrestrial springs in Slovakia. Additionally, the study emphasizes that sulfur springs in Slovakia support unique, poorly characterized bacterial communities of sulfur-oxidizing bacteria.

Profiling of adrenal corticosteroids in blood and local tissues of mice during chronic stress.

Stress increases plasma concentrations of corticosteroids, however, their tissue levels are unclear. Using a repeated social defeat paradigm, we examined the impact of chronic stress on tissue levels of corticosterone (CORT), progesterone (PROG), 11-deoxycorticosterone (11DOC) and 11-dehydrocorticosterone (11DHC) and on gut microbiota, which may reshape the stress response. Male BALB/c mice, liquid chromatography-tandem mass spectrometry and 16S RNA gene sequencing were used to screen steroid levels and fecal microbiome, respectively. Stress induced greater increase of CORT in the brain, liver, and kidney than in the colon and lymphoid organs, whereas 11DHC was the highest in the colon, liver and kidney and much lower in the brain and lymphoid organs. The CORT/11DHC ratio in plasma was similar to the brain but much lower in other organs. Stress also altered tissue levels of PROG and 11DOC and the PROG/11DOC ratio was much higher in lymphoid organs that in plasma and other organs. Stress impacted the ß- but not the α-diversity of the gut microbiota and LEfSe analysis revealed several biomarkers associated with stress treatment. Our data indicate that social defeat stress modulates gut microbiota diversity and induces tissue-dependent changes in local levels of corticosteroids, which often do not reflect their systemic levels.

Host Species Affects Bacterial Evenness, but Not Diversity: Comparison of Fecal Bacteria of Cows and Goats Offered the Same Diet.

The aim of this study was to compare the diversity and composition of fecal bacteria in goats and cows offered the same diet and to evaluate the influence of animal species on the gut microbiome. A total of 17 female goats (Blond Adamellan) and 16 female cows (Brown Swiss) kept on an organic farm were fed pasture and hay. Bacterial structure in feces was examined by high-throughput sequencing using the V4-V5 region of the 16S rRNA gene. The Alpha diversity measurements of the bacterial community showed no statistical differences in species richness and diversity between the two groups of ruminants. However, the Pielou evenness index revealed a significant difference and showed higher species evenness in cows compared to goats. Beta diversity measurements showed statistical dissimilarities and significant clustering of bacterial composition between goats and cows. Firmicutes were the dominant phylum in both goats and cows, followed by Bacteroidetes, Proteobacteria, and Spirochaetes. Linear discriminant analysis with effect size (LEfSe) showed a total of 36 significantly different taxa between goats and cows. Notably, the relative abundance of Ruminococcaceae UCG-005, Christensenellaceae R-7 group, Ruminococcaceae UCG-010, Ruminococcaceae UCG-009, Ruminococcaceae UCG-013, Ruminococcaceae UCG-014, Ruminococcus 1, Ruminococcaceae UCG-002, Lachnospiraceae NK4A136 group, Treponema 2, Lachnospiraceae AC2044 group, and Bacillus was higher in goats compared to cows. In contrast, the relative abundance of Turicibacter, Solibacillus, Alloprevotella, Prevotellaceae UCG-001, Negativibacillus, Lachnospiraceae UCG-006, and Eubacterium hallii group was higher in cows compared with goats. Our results suggest that diet shapes the bacterial community in feces, but the host species has a significant impact on community structure, as reflected primarily in the relative abundance of certain taxa.

Dysbiosis of skin microbiome and gut microbiome in melanoma progression.

BACKGROUND: The microbiome alterations are associated with cancer growth and may influence the immune system and response to therapy. Particularly, the gut microbiome has been recently shown to modulate response to melanoma immunotherapy. However, the role of the skin microbiome has not been well explored in the skin tumour microenvironment and the link between the gut microbiome and skin microbiome has not been investigated in melanoma progression. Therefore, the aim of the present study was to examine associations between dysbiosis in the skin and gut microbiome and the melanoma growth using MeLiM porcine model of melanoma progression and spontaneous regression. RESULTS: Parallel analysis of cutaneous microbiota and faecal microbiota of the same individuals was performed in 8 to 12 weeks old MeLiM piglets. The bacterial composition of samples was analysed by high throughput sequencing of the V4-V5 region of the 16S rRNA gene. A significant difference in microbiome diversity and richness between melanoma tissue and healthy skin and between the faecal microbiome of MeLiM piglets and control piglets were observed. Both Principal Coordinate Analysis and Non-metric multidimensional scaling revealed dissimilarities between different bacterial communities. Linear discriminant analysis effect size at the genus level determined different potential biomarkers in multiple bacterial communities. Lactobacillus, Clostridium sensu stricto 1 and Corynebacterium 1 were the most discriminately higher genera in the healthy skin microbiome, while Fusobacterium, Trueperella, Staphylococcus, Streptococcus and Bacteroides were discriminately abundant in melanoma tissue microbiome. Bacteroides, Fusobacterium and Escherichia-Shigella were associated with the faecal microbiota of MeLiM piglets. Potential functional pathways analysis based on the KEGG database indicated significant differences in the predicted profile metabolisms between the healthy skin microbiome and melanoma tissue microbiome. The faecal microbiome of MeLiM piglets was enriched by genes related to membrane transports pathways allowing for the increase of intestinal permeability and alteration of the intestinal mucosal barrier. CONCLUSION: The associations between melanoma progression and dysbiosis in the skin microbiome as well as dysbiosis in the gut microbiome were identified. Results provide promising information for further studies on the local skin and gut microbiome involvement in melanoma progression and may support the development of new therapeutic approaches.

Bifidobacterium canis sp. nov., a novel member of the Bifidobacterium pseudolongum phylogenetic group isolated from faeces of a dog (Canis lupus f. familiaris).

A fructose-6-phosphate phosphoketolase-positive strain (GSD1FST) was isolated from a faecal sample of a 3 weeks old German Shepherd dog. The closest related taxa to isolate GSD1FST based on results from the EZBioCloud database were Bifidobacterium animalis subsp. animalis ATCC 25527T, Bifidobacterium animalis subsp. lactis DSM 10140T and Bifidobacterium anseris LMG 30189T, belonging to the Bifidobacterium pseudolongum phylogenetic group. The resulting 16S rRNA gene identities (compared length of 1454 nucleotides) towards these taxa were 97.30, 97.23 and 97.09 %, respectively. The pairwise similarities of strain GSD1FST using argS, atpA, fusA, hsp60, pyrG, rpsC, thrS and xfp gene fragments to all valid representatives of the B. pseudolongum phylogenetic group were in the concatenated range of 83.08-88.34 %. Phylogenomic analysis based on whole-genome methods such as average nucleotide identity revealed that bifidobacterial strain GSD1FST exhibits close phylogenetic relatedness (88.17 %) to Bifidobacetrium cuniculi LMG 10738T. Genotypic characteristics and phylogenetic analyses based on nine molecular markers, as well as genomic and comparative phenotypic analyses, clearly proved that the evaluated strain should be considered as representing a novel species within the B. pseudolongum phylogenetic group named as Bifidobacterium canis sp. nov. (GSD1FST=DSM 105923T=LMG 30345T=CCM 8806T).

Glutamine synthetase type I (glnAI) represents a rewarding molecular marker in the classification of bifidobacteria and related genera.

The family Bifidobacteriaceae constitutes an important phylogenetic group that particularly includes bifidobacterial taxa demonstrating proven or debated positive effects on host health. The increasingly widespread application of probiotic cultures in the twenty-first century requires detailed classification to the level of particular strains. This study aimed to apply the glutamine synthetase class I (glnAI) gene region (717 bp representing approximately 50% of the entire gene sequence) using specific PCR primers for the classification, typing, and phylogenetic analysis of bifidobacteria and closely related scardovial genera. In the family Bifidobacteriaceae, this is the first report on the use of this gene for such purposes. To achieve high-value results, almost all valid Bifidobacteriaceae type strains (75) and 15 strains isolated from various environments were evaluated. The threshold value of the glnAI gene identity among Bifidobacterium species (86.9%) was comparable to that of other phylogenetic/identification markers proposed for bifidobacteria and was much lower compared to the 16S rRNA gene. Further statistical and phylogenetic analyses suggest that the glnAI gene can be applied as a novel genetic marker in the classification, genotyping, and phylogenetic analysis of isolates belonging to the family Bifidobacteriaceae.

Genetic marker-based multi-locus sequence analysis for classification, genotyping, and phylogenetics of the family Bifidobacteriaceae as an alternative approach to phylogenomics.

Bifidobacteria are widely known for their probiotic potential; however, little is known regarding the ecological significance and potential probiotic effects of the phylogenetically related 'scardovial' genera (Aeriscardovia, Alloscardovia, Bombiscardovia, Galliscardovia, Neoscardovia, Parascardovia, Pseudoscardovia and Scardovia) and Gardnerella classified with bifidobacteria within the Bifidobacteriaceae family. Accurate classification and genotyping of bacteria using certain housekeeping genes is possible, whilst current phylogenomic analyses allow for extremely precise classification. Studies of applicable genetic markers may provide results comparable to those obtained from phylogenomic analyses of the family Bifidobacteriaceae. Segments of the glyS (624 nucleotides), pheS (555 nucleotides), rpsA (630 nucleotides), and rpsB (432 nucleotides) genes and their concatenated sequence were explored. The mean glyS, pheS, rpsB and rpsA gene sequence similarities calculated for Bifidobacterium taxa were 84.8, 85.2, 90.2 and 86.8%, respectively. Interestingly, the average value of the Average Nucleotide Identity among 67 type strains of the family Bifidobacteriaceae (84.70%) calculated based on values published recently was in agreement with the average pairwise similarity (84.6%) among 75 type strains of Bifidobacteriaceae family computed in this study using the concatenated sequences of four gene fragments. Similar to phylogenomic analyses, several gene sequence and phylogenetic analyses revealed that concatenated gene regions allow for classification of Bifidobacteriaceae strains into particular phylogenetic clusters and groups. Phylogeny reconstructed from the concatenated sequences assisted in defining two novel phylogenetic groups, the Bifidobacterium psychraerophilum group consisting of B. psychraerophilum, Bifidobacterium crudilactis and Bifidobacterium aquikefiri species and the Bifidobacterium bombi group consisting of B. bombi, Bifidobacterium bohemicum and Bifidobacterium commune.

Effect of Selected Stilbenoids on Human Fecal Microbiota.

Dietary phenolics or polyphenols are mostly metabolized by the human gut microbiota. These metabolites appear to confer the beneficial health effects attributed to phenolics. Microbial composition affects the type of metabolites produced. Reciprocally, phenolics modulate microbial composition. Understanding this relationship could be used to positively impact health by phenolic supplementation and thus create favorable colonic conditions. This study explored the effect of six stilbenoids (batatasin III, oxyresveratrol, piceatannol, pinostilbene, resveratrol, thunalbene) on the gut microbiota composition. Stilbenoids were anaerobically fermented with fecal bacteria from four donors, samples were collected at 0 and 24 h, and effects on the microbiota were assessed by 16S rRNA gene sequencing. Statistical tests identified affected microbes at three taxonomic levels. Observed microbial composition modulation by stilbenoids included a decrease in the Firmicutes to Bacteroidetes ratio, a decrease in the relative abundance of strains from the genus Clostridium, and effects on the family Lachnospiraceae. A frequently observed effect was a further decrease of the relative abundance when compared to the control. An opposite effect to the control was observed for Faecalibacterium prausnitzii, whose relative abundance increased. Observed effects were more frequently attributed to resveratrol and piceatannol, followed by thunalbene and batatasin III.

Melanoma-related changes in skin microbiome.

Melanoma is the least common form of skin tumor, but it is potentially the most dangerous and responsible for the majority of skin cancer deaths. We suggest that the skin microbiome might be changed during the progression of melanoma. The aim of this study is to compare the composition of the skin microbiota between different locations (skin and melanoma) of a MeLiM (Melanoma-bearing Libechov Minipig) pig model (exophytic melanoma). Ninety samples were used for PCR-DGGE analysis with primers specifically targeting the V3 region of the 16S rRNA gene. The profiles were used for cluster analysis by UPGMA and principal coordinate analysis PCoA and also to calculate the diversity index (Simpson index of diversity). By comparing the obtained results, we found that both bacterial composition and diversity were significantly different between the skin and melanoma microbiomes. The abundances of Fusobacterium and Trueperella genera were significantly increased in melanoma samples, suggesting a strong relationship between melanoma development and skin microbiome changes.

Fragment of the aspartyl-tRNA synthetase applicable as a shared classification and phylogenetic marker in particular representatives of the order Lactobacillales.

The order Lactobacillales represents a morphologically, metabolically, and physiologically diverse group of bacteria. Lactic acid bacteria represent the core of this phylogenetic group. They are a part of epiphytic microflora, fermented dairy, meat, fruit and vegetable products, and the digestive tract of humans and animals. Despite the fact that these bacteria form a phenotypically and genotypically heterogeneous group, their phylogenetic relationship enables to propose a common genetic marker usable in classification, typing, and phylogeny. By creation of consensus sequence based on available genomic sequences of some representatives of order Lactobacillales, a specific primer-pair binding variable region of aspS gene (length of 615 nts) encoding the aspartyl-tRNA synthetase was designed. This gene has not yet been used in classification and phylogeny of the order Lactobacillales, although it meets the requirements of molecular markers (distribution and single copy in bacterial genomes, functional constancy and genetic stability, sequence variability among taxonomic units, irreplaceable role in proteosynthesis). Primers were applied on 54 type and wild Lactobacillales strains. Obtained sequences allowed to provide alignments for purpose of phylogenetic tree reconstructions that uncovered particular phylogenetic clusters of vagococci/enterococci, obligately homofermentative and heterofermentative lactobacilli. Although a relatively short fragment of the aspS gene (approximately 33% of the complete gene sequence) was evaluated, much higher sequence variability (61.8% of pairwise identity) among strains examined compared with 16S rRNA gene (90.7%, length of 1318 nt) provides a relatively simple and effective tool for classification and typing of selected representatives of the order Lactobacillales.

The threonine-tRNA ligase gene region is applicable in classification, typing, and phylogenetic analysis of bifidobacteria.

In the modern era, molecular genetic techniques are crucial in ecological studies, as well as in the classification, typing, and phylogenetic analysis of prokaryotes. These techniques are mainly aimed at whole genome comparisons and PCR-derived experiments, including amplifying the 16S rRNA and other various housekeeping genes used in taxonomy, as well as MLST (multilocus sequence typing) and MLSA (multilocus sequence analysis) of different taxonomic bacterial groups. The gene encoding threonine-tRNA ligase (thrS) is a gene potentially applicable as an identification and phylogenetic marker in bacteria. It is widely distributed in bacterial genomes and is subject to evolutionary selection pressure due to its important function in protein synthesis. In this study, specific primers were used to amplify a thrS gene fragment (~740 bp) in 36 type and 30 wild strains classified under family Bifidobacteriaceae. The full-length gene has not yet been considered as a possible identification, classification, and phylogenetic marker in bifidobacteria. The thrS sequences revealed higher sequence variability (82.7% of pairwise identities) among members of the family than that shown by 16S rRNA gene sequences (96.0%). Although discrepancies were found between the thrS-derived and previously reported whole genome phylogenetic analyses, the main phylogenetic groups of bifidobacteria were properly assigned. Most wild strains of bifidobacteria were better differentiated based on their thrS sequences than on their 16S rRNA gene identities. Phylogenetic confidence of the evaluated gene with respect to other alternative genetic markers widely used in taxonomy of bifidobacteria (fusA, GroELhsp60, pyrG, and rplB genes) was confirmed using the localized incongruence difference - Templeton analysis.

Gene encoding the CTP synthetase as an appropriate molecular tool for identification and phylogenetic study of the family Bifidobacteriaceae.

An alternative molecular marker with respect to the 16S rRNA gene demonstrating better identification and phylogenetic parameters has not been designed for the whole Bifidobacteriaceae family, which includes the genus Bifidobacterium and scardovial genera. Therefore, the aim of the study was to find such a gene in available genomic sequences, suggest appropriate means and conditions for asmplification and sequencing of the desired region of the selected gene in various strains of the bacterial family and verify the importance in classification and phylogeny. Specific primers flanking the variable region (~800 pb) within the pyrG gene encoding the CTP synthetase were designed by means of gene sequences retrieved from the genomes of strains belonging to the family Bifidobacteriaceae. The functionality and specificity of the primers were subsequently tested on the wild (7) and type strains of bifidobacteria (36) and scardovia (7). Comparative and phylogenetic studies based on obtained sequences revealed actual significance in classification and phylogeny of the Bifidobacteriaceae family. Gene statistics (percentages of mean sequence similarities and identical sites, mean number of nucleotide differences, P- and K-distances) and phylogenetic analyses (congruence between tree topologies, percentages of bootstrap values >50 and 70%) indicate that the pyrG gene represents an alternative identification and phylogenetic marker exhibiting higher discriminatory power among strains, (sub)species, and genera than the 16S rRNA gene. Sequences of the particular gene fragment, simply achieved through specific primers, enable more precisely to classify and evaluate phylogeny of the family Bifidobacteriaceae including, with some exceptions, health-promoting probiotic bacteria.

Alloscardovia venturai sp. nov., a fructose 6-phosphate phosphoketolase-positive species isolated from the oral cavity of a guinea-pig (Cavia aperea f. porcellus).

A slightly irregular, short rod-shaped bacterial strain, MOZIV/2T, showing activity of fructose 6-phosphate phosphoketolase was isolated from the oral cavity of a home-bred guinea-pig. Based on comparative 16S rRNA gene sequence analyses, its closest relatives were Alloscardovia omnicolens DSM 21503T and Alloscardovia criceti DSM 17774T with 96.0 and 95.6 % pairwise similarities, respectively. Completeness of the compared sequences was 97.3 and 96.9 %, respectively. Growth was found only under anaerobic conditions. Activities of α- and ß-gluco(galacto)sidases were detected in strain MOZIV/2T, which is characteristic for almost all members of the family Bifidobacteriaceae. Sequencing of other molecular markers (fusA, gyrB and xfp) revealed low gene sequence similarities to A. omnicolens DSM 21503T ranging from 72.7 to 87.5 %. Strain MOZIV/2T differed from other species within the genus Alloscardovia by the presence of C18 : 1ω9t. In addition, much higher proportions of C8 : 0, C11 : 0, C12 : 0, C14 : 1, C16 : 1 and C17 : 0 fatty acids were found in cells of strain MOZIV/2T. The peptidoglycan structure was of type A4α [l-Lys(l-Orn)-d-Asp], which is consistent with its classification within the genus Alloscardovia. The DNA G+C content (45.8 mol%) was lower than those found in other alloscardovia. Phylogenetic studies and evaluation of phenotypic characteristics including the results of biochemical, physiological and chemotaxonomic analyses confirmed the novel species status for strain MOZIV/2T, for which the name Alloscardovia venturai sp. nov. is proposed. The type strain is MOZIV/2T (=DSM 100237T=CCM 8604T=LMG 28781T).

Lactobacillus caviae sp. nov., an obligately heterofermentative bacterium isolated from the oral cavity of a guinea pig (Cavia aperea f. porcellus).

A Gram-stain-positive, facultatively anaerobic, and catalase- and oxidase-negative bacterial strain designated MOZM2T, having 98.4 % 16S rRNA gene sequence identity with Lactobacillus reuteri DSM 20016T, was isolated from a swab of the oral cavity of a home-bred guinea pig. Comparative analyses based on the hsp60, pheS and tuf genes confirmed L. reuteri as its closest relative species, with calculated sequence similarities of 92.8, 88.8 and 96.9 %, respectively. DNA-DNA hybridisation revealed a 42 % degree of genetic similarity between the novel strain and L. reuteri DSM 20016T. Strain MOZM2T degrades carbohydrates via the 6-phosphogluconate/phosphoketolase pathway, evidenced by its production of gas from glucose and the end products of hexose catabolism. Comparative analysis of the cellular fatty acid profiles determined significant differences between MOZM2T and L. reuteri DSM 20016T in their proportions of C8 : 0, C14 : 1, C17 : 0, C18 : 2ω6t and C20 : 0 fatty acids. Results of genotypic analyses also demonstrated differences between these two strains. They also differed in DNA G+C content, and some biochemical and physiological characteristics. We therefore believe that the examined bacterial isolate should be considered as a new taxon within the group of obligately heterofermentative lactobacilli. The species name Lactobacillus caviae sp. nov. is proposed, of which the type strain is MOZM2T (=CCM 8609T=DSM 100239T=LMG 28780T).

Classification of Culturable Bifidobacterial Population from Colonic Samples of Wild Pigs (Sus scrofa) Based on Three Molecular Genetic Methods.

Occurrence of bifidobacteria, known as health-promoting probiotic microorganisms, in the digestive tract of wild pigs (Sus scrofa) has not been examined yet. One hundred forty-nine fructose-6-phosphate phosphoketolase positive bacterial strains were isolated from colonic content of twenty-two individuals of wild pigs originated from four localities in the Czechia. Based on PCR-DGGE technique targeting the variable V3 region of the 16S rRNA genes, strains were initially differentiated into four groups represented by: (i) probably a new Bifidobacterium species (89 strains), (ii) B. boum/B. thermophilum/B. thermacidophilum subsp. porcinum/B. thermacidophilum subsp. thermacidophilum (sub)species (49 strains), (iii) Pseudoscardovia suis (7 strains), and (iv) B. pseudolongum subsp. globosum/B. pseudolongum subsp. pseudolongum (4 strains), respectively. Given the fact that DGGE technique did not allow to differentiate the representatives of thermophilic bifidobacteria and B. pseudolongum subspecies, strains were further classified by the 16S rRNA and thrS gene sequences. Primers targeting the variable regions of the latter gene were designed to be applicable in identification and phylogeny of Bifidobacteriaceae family. The 16S rRNA-derived phylogenetic study classified members of the first group into five subgroups in a separated cluster of thermophilic bifidobacteria. Comparable results were obtained by the thrS-derived phylogenetic analysis. Remarkably, variability among thrS sequences was higher compared with 16S rRNA gene sequences. Overall, molecular genetic techniques application allowed to identify a new Bifidobacterium phylotype which is predominant in the digestive tract of examined wild pigs.